RUNHETC-2000-17 
hepth / nnnmmyy 



Finite Temperature Expectation Values 
of Local Fields in the sinh-Gordon model 



Sergei Lukyanov 

Department of Physics and Astronomy, Rutgers University 
Piscataway, NJ 08855-0849, USA 
and 

L.D. Landau Institute for Theoretical Physics 
Kosygina 2, Moscow, Russia 



Abstract 

Sklyanin's method of separation of variables is employed in a calculation of finite tem- 
perature expectation values. An essential element of the approach is Baxter's Q-function. 
We propose its explicit form corresponding to the ground state of the sinh-Gordon theory 
With the method of separation of variables we calculate the finite temperature expectation 
values of the exponential fields to one-loop order of the semi-classical expansion. 
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1. Introduction 

One-point functions have numerous applications in Statistical Mechanics and Con- 
densed Matter Physics They determine "generalized susceptibilities" i.e. the linear 
response of a system to external fields. In a path integral formulation the one-point func- 
tion of a local field O is represented by a Euclidean path integral of the form 

-A 



(O) = Z~ L J T>[cp] O e~ A . (1.1) 

Recently some progress has been achieved in the calculation of one-point functions in 
integrable Quantum Field Theory (QFT) defined on a two dimensional Euclidean plane 
@,(§. In this case, the integral ( p..l| ) can also be viewed as a Vacuum Expectation Value 
(VEV) in 1 + 1-dimensional QFT associated with the action A. For many applications, 
especially in Condensed Matter Physics, it is important to generalize the results of Refs.|||] 
to the case of Euclidean path integrals defined on an infinite cylinder. In the Matsubara 
imaginary time formalism such path integrals are interpreted as thermal averages 

where H is the Hamiltonian of the corresponding QFT and the temperature coincides with 
the inverse circumference of the cylinder. 

The path integral defined on a cylinder also allows another interpretation. It is an 
expectation value for the ground state | vac)n of the 1 + 1-dimensional theory in the finite 
geometry where the spatial coordinate is compactified on a circle. Hence, the VEVs contain 
important information about Renormalization Group flow controlled by the parameter R. 

An exact calculation of the finite volume (finite R) VEVs is a challenge even in in- 
tegrable QFT. Recent progress made in papers should be mentioned here. In || A. 
Leclair and G. Mussardo proposed an integral representation which makes it possible to 
generate a low-temperature (R — > oo) expansion for the VEVs in terms of infinite vol- 
ume form-factors of local fields and some thermodynamical data. Their conjecture works 
for theories with trivial S'-matrices such as Ising and Free Dirac Fermion models 0, but 
its validity remains questionable for models with non-trivial scattering amplitudes (see 
e.g. H). Another line of research was proposed in the work ||. F. Smirnov applied the 
method of separation of variables PJTO|] to the semi-classical study of finite volume matrix 
elements in the quantum KdV equation. The model does not constitute a relativistic field 
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theory. Nevertheless, it is of prime importance for the sinh-Gordon QFT since both of the 
equations are in the same integrable hierarchy. 

In this paper we will implement the method of separation of variables in the case of 
the quantum sinh-Gordon theory. The problem is defined by a Euclidean action, 

AhG = J dxi J ^{J-(^) 2 + 2,c OS h(M}, (I- 3 ) 

where (p is a scalar field with periodic boundary condition along the ^-coordinate. We 
are focusing on the VEVs of the exponential fields, 

O = e atp . (1.4) 

For our purposes it will be useful to rewrite (1.2[) in the form, 



(e a n R = Z- 1 jv[x]^l[x]e ax ■ (1-5) 

Here VPo[x] is an integral taken over field configurations on the half-cylinder, x\ < 0, 
satisfying the boundary condition, 

i.e. it is a wave functional corresponding to the ground state \ vo,c)r. The method of 
separation of variables |J],[T0|,|5[] allows one to introduce a change of integration variables in 

(ED, 

xm -> {7fc}r=-oc' 

from the function x{ x 2) to the infinite discrete set of 7^. A notable advantage of the new 
variables is that the wave functional in the "7-representation" has a factorizable form, 

*o[{7a] ~n s ^ • 

k 

Notice that the integration measure ^[{7/c}] does not factorize in the variables 7^. At 
this moment, we are not able to elaborate on all steps of the changing of variables on 
a rigorous basis. Therefore, we suggest the deduced integral representation for (e a{p )n 
as a conjecture rather than a well established result. To test the validity of this integral 
representation, we carry out a semi-classical expansion of the VEV. 



2 



More explicitly, the parameter b 2 in the action (|1.3| ) can readily be identified with 
the Planck constant. Then, for finite a = ab and b 2 — > 0, the functional integral ( |1 . 1| ) 
is dominated by a non-trivial saddle-point configuration and admits the semi-classical 
expansion, 

(e a «') R = e-$ D (l + 0(6 2 )) . (1.6) 

Here S is a Euclidean action on the cylinder evaluated in the saddle-point configuration 
and the pre-exponential factor D is the result of evaluating the functional integral ( |1 . 1| ) 
in the Gaussian approximation around the classical solution. With the proposed integral 
representation we calculate the functions S and D and find complete agreement with 
the expected high- and low-temperature behavior of the VEVs. In particular, our result 
matches well with the Leclair-Mussardo conjecture. 



2. Integral representation for VEVs 

2.1. Flaschka- McLaughlin variables 



In the paper |Tl| H. Flaschka and D. McLaughlin found remarkable canonically con- 
jugate variables in the phase spaces of the classical Toda chain and KdV equations. Their 
approach can be straightforwardly adapted to the classical sinh-Gordon equation. Here 
we give a brief review of the Flaschka-McLaughlin variables for this model. For more 
information and proofs, the reader is referred to Refe.[H,|I2|. 



The sinh-Gordon equation admits a zero curvature formulation: There exists a 
si ( 2, R)- valued connection 1-form, depending on an auxiliary parameter A, such that the 
condition of vanishing curvature is equivalent to the equation of motion. Dealing with the 
theory on cylinder, one can integrate this 1-form along some cycle, say, 

x 1 = 0, 0<x 2 <R, (2.1) 

and obtain the so-called monodromy matrix, 

M{\)=(*&) ^[y] e SL(2,M) (3mA = 0). (2.2) 



A" 1 C(A) D(A) 

This matrix satisfies important analytical conditions which are readily obtained from an 
explicit form of the connection. In particular, the matrix elements in ( |2.2| ) are real analyt- 
ical functions of the variable A 2 with two essential singularities at the points A 2 = 0, oo. 
Zeroes of B(A), 

A 2 : B(A fc ) = , (2.3) 
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are of prime importance in the construction. It is possible to show that all zeroes ( |2.3| ) 
are simple, real, positive and accumulate towards the essential singularities. Thus we can 
order them, 

«- ...X 2 _ N < X- N+1 . . . < Aq < . < x n ■ ■ ■ -» +°o 

and define two infinite sets: 

{7fc}^-oo ■■ Ik = logA^, 

oo" ( 2 - 4 ) 
WJfcl-oo : ^ = 4 lo Si A ( A fe)| ■ 

The mapping of the canonical Poisson data, (p and d Xl (p , to the variables (|2.4|) is found to 
be a canonical transformation, i.e. 

{ TTfc , 7m } = 4m , { ^ k , } = { Ik , 1m } = ■ (2.5) 

Hence ( |2.4j ) are canonically conjugate variables in the phase space of the sinh-Gordon 
model. At the same time, we can treat {lk}'^ 3 __ as coordinates in the corresponding 
configuration space. 

2.2. 7 -representation 

The Flaschka-McLaughlin variables proved to be useful in quantum theory as was 
demonstrated in the seminal work || on the example of the Toda chain equation (see also 
||T^,[14|)- We refer to a quantization in these variables as a quantization in 7-representation. 
Recently the 7-representation was employed to quantize "real" KdV theory ||0. In fact, 
Smirnov presented heuristic, but convincing, mo del- independent arguments which can also 
be applied to the sinh-Gordon equation. Following these arguments we introduce the 
integral, 

N /.+00 7 N 

Z N (R,a)= U II sinh( 7fc - 7m ) II Q 2 h*} e 2 -^ (ab+k) ■ (2-6) 

fc=-iV J 00 N>k>m>-N k=-N 

1 The monodromy matrix in the "real" KdV model has the form ( |2.2| ), whereas the monodromy 
matrix of the "imaginary" equation belongs to the group SU(2) (9mA = 0). The imaginary KdV 
model is related to the sine-Gordon theory and (perturbed) CFT with the central charge c < 1 
(see e.g. [p!5[). A sensible 7-representation for the imaginary equation has, to our knowledge, not 
been found. 
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We shall also use slightly different form of In '■ With the identity. 

n 2 S in h (2i^). Det |exp(^ 



b 2 

N>k>m>-N 



-N<j,k<N 



the integral ( |2.6| ) can be rewritten in the form 

1 N f 



N 

JI 2 sinh( 7fc - 7m ) sinh ( ) J] ^\ ^ 

N>k>m>-N " k=-N 



(2.7) 



The function Q[y] appearing in Eqs. (|2.6|) , (|2.7| ) is the so-called Baxter's Q-function. It is 
a non-singular function for real 7 with leading asymptotic behavior 

Q[ 7 ] ~ e" 2Co cosh (^T) as 7 -> 00 . (2.8) 

Here and below we use the notation, 

q = b + b~ l . 



The positive constant Cq in ( |2.8|) reads explicitly, 

mi? 

^0 



4 sin(^) 



and m is a mass of the sinh-Gordon particle. Therefore, In Q2.7 ) is a convergent integral 



for any finite N. Notice that the configuration space of the model under consideration 
is an infinite-dimensional space. In writing the 2N + 1-fold integral, we truncate it to 
2N + 1-dimensional space with coordinates {jk} k= _ N - ^ s well as in the KdV theory ||, 
we can treat 

N 

**[{7*}] = n 

fc=-JV 

as a wave functional in the 7- representation. For iV — > 00, it corresponds to the ground 
state I vac)n of the sinh-Gordon theory with periodic boundary conditions. Furthermore, 
the double product in ( |2.7| ) is an integration measure and we shall denote it Pjv[{7}]- 
It was calculated in the semi-classical approximation in 0. The semi-classical analysis 
suggests also that the product 



o N = n 



N 

2 ~<k a 

e *> 



k=-N 
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represents (in the limit N — > oo) the exponential field e a<p located at the point (-R/2, R/2) 
on the cylinder!. Therefore, the integral ( |2.7|) has the form of a quantum mechanical 
diagonal matrix element, 

Tn = Jv n [ {7} ] ^ N O n q N . (2.9) 

We shall consider the Vacuum Expectation Values only. In this case, Q is an eigenvalue of 
the Baxter Q-operator corresponding to the ground state. 

2.3. Baxter's Q -function in sink- Gordon theory 

We now turn to the most delicate point of our construction: an explicit form of the 
function Q. To the best of our knowledge a rigorous derivation of Q is not currently 
available. Here we formulate a conjecture for the sinh-Gordon Q-function based on the 
following heuristic arguments. First, we note that the substitution b —>■ i (3 transforms *A s hG 
( |1.3|) to the action of the sine-Gordon model. Naively we could try to obtain Q by means 
of analytical continuation from — 1 < b 2 < 0. In this coupling constant domain, the Baxter 
Q-operator is relatively well studied (|T5fl . Unfortunately, the sine-Gordon Q- function has 
an essential singularity at b 2 = the analytical structure of which is unknown. This makes 
the analytical continuation to the domain of the sinh-Gordon model a highly questionable 
procedure. One can guess an explicit form for Q by examining its asymptotic behavior. 
The Q-function in the sine-Gordon model admits the following asymptotic expansion for 
7 -> +00 pflll, 

00 00 
logQ sinG ~ -C e 9 + Y,Cn hn-ie-^-W + £ G n e 2 ^^ (-1 < b 2 < 0) . 

ra=l n=l 

Here we introduce a new notation, 

9 = jbq. (2.10) 

The leading term of this expansion has already appeared in our consideration (see ( |2.8| )) 
while I271-1 and G n are vacuum eigenvalues of the so-called local and dual unlocal Integrals 
of Motion (IM) respectively. The constants C n depend on the normalization of the local IM. 

2 The position of the insertion is determined by choosing of the integration contour for the 
monodromy matrix. Here we assume that the contour is given by ( |2.1| ). 

3 In this work we use a convention for the Q-function which differs from the one of [ fL5| by an 
overall shift of the argument. 
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In Appendix A (see Eq.(A.2)) we present their form for the normalization adopted in [15 



A similar asymptotic form holds for 7 — ► —00. The eigenvalues l2n-i are regular functions 
of b 2 , and can be continued to the domain of the sinh-Gordon model without problems. 
Contrary to l2n-i, the eigenvalues of the dual unlocal IM, G n , are highly singular functions 
at b 2 = 0. One can expect that the appearance of these IM is a consequence of the existence 
of the soliton sector of the sine-Gordon QFT. This sector is absent for b 2 > 0. All these 
observations suggest the following large 7 asymptotic behavior in the sinh-Gordon model, 

00 

log Q ~ -Co e e - C n I 2n _i e -C 2 »-D* (b 2 > 0) . (2.11) 

71=1 

In Appendix A we give numerical evidence that the values of local IM can be expressed in 
terms of a solution of the Thermodynamic Bethe Ansatz (TBA) equation: 

C n I 2n _! = Co S nl + (-1T J°° ® log ( 1 + e-^ ) . (2.12) 



Here the function e(9) solves the TBA equation [ [[6|jl7jjl8| 

/oo mi . . 

— <S>(e-6') logfl + e-^ 9 ') I =0 , (2.13) 
-00 2tt V J 

with the kernel 

4 sin (2^) cosh(^) 

$(6) = t^4t • (2-14) 

1 ' cosh(2^) -cos(^) 1 ' 

The series (|2.11| ) is an asymptotic expansion. In fact, it is a divergent geometrical series 
which can easily be summed up and one can guess an explicit form of Q: 

log m = cosm + r * Mj ;«-;''') . ( , 16) 

2 sm (y) J-oo 27T cosh(6' - 0') 

We leave an examination of the properties of this function for future publications. 

One more aspect of the 7-representation deserves a comment. The sinh-Gordon model 
manifests an important non-perturbative symmetry. The couplings b and b~ l correspond 



to physically indistinguishable theories. Q in (|2.15|) is a self-dual function in a sense that 
it is invariant under the substitution b —>■ Furthermore, it is easy to see that 

2jvL = 2jvL_i ■ 



This supports our choice of the measure in (|2.7|) . Strictly speaking, the measure was 
obtained in |5[ in the semi-classical approximation. The exact invariance of the semi- 
classical measure suggests its applicability for an arbitrary value of the coupling constant 
b 2 . 
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2.4- Large N limit 

As was noted above, in writing the 2 N + 1-fold integral ( p.7|) , we truncate the con- 
figuration space of the theory to 2N + 1-dimensional one. In fact, the truncation amounts 
to an ultraviolet regularization with a momentum cutoff given by 



A;v = 



2nN 
R 



Now we let iV — > oo. From (|2.9|) , it is clear that the ratio 



1n(R, a) 



Zn(R,0) 



(2.16) 



represents the VEV {e a(p )n in the large N limit. More explicitly, dimensional analysis 
suggests that 

Zn(R, a) oc 

V m 



thus the correct relation has the form 



(e a *) R = K a lim^ 00 (— ) 1 N (R, 



(2.17) 



Here n a is an arbitrary i?-independent constant. To eliminate this ambiguity one has 
to impose some normalization condition on the fields. For example, the so-called con- 
formal normalization stipulates that the exponential fields with sufficiently small \a\ are 
normalized in accordance with the short distance behavior of the two-point function 



(e^(x)e-^(y))^\x-y 



4a" 



as \x — y \ . 



The key result of Refs. |3|,|J is a calculation of the limit, 



hm (e at ?) R = g a 

R — »oo 



(2.18) 



with this normalization. Explicitly, 



Ga 



m . 



r(i)r(i + £) 



2q, 



exp 



00 dt 

T 



-2a z 



sinh 2 (2a6t) 



2 sinh(6 2 t) sinh(t) cosh.(qbt) 



+ 2a 2 e~ 2t 



(2.19) 



Once we adopt the conformal normalization for the exponential fields, the constant K a is 
uniquely determined by the condition Q2.18Q . In particular, the semi-classical consideration 
(see below) leads to the relation, 

Ka = G a (l + 0(6 2 )) . 



3. Semi-classical expansion 

In this section we will study the VEVs in the semi-classical approximation. The VEV 
( P-17| ) can be represented by the Euclidean path integral on a cylinder with A and 
O given by (|L3l) , For fixed a, 

a = ab (3-1) 

and b 2 — > 0, the path integral is dominated by a saddle-point configuration (j) = btp and 
the VEV has the form ( |1.6| ), where S coincides with the regularized Euclidean classical 
action on the cylinder: 

d 2 x { (drrd)) 2 n . , o /6\ 1 f ds 



S 



lim 



8tt 



(d a (j>) o 2 

Ym sinh , 

2 V2 



+ a 



\x-y\=e 27r 



<-2a loge 



\x-y\>e 

Here the function <p is a solution of the classical equation of motion 

d 2 4> = m 2 sinh((/>) , 



such that, 



and 



4a log \x — y\ + 0(1) as \x — y\ 







(3.2) 

(3.3) 
(3.4) 



y\ 



oo 



(3.5) 



4> — > as \x 

4>(xi,x 2 + R) = 4>(xi,x 2 ) ■ 
The field configuration (p develops a singularity at the point y where the exponential field 
is inserted. Therefore, in the definition ( |3.2|) we cut the small disc of radius e around 
this point and add the boundary term to the action to ensure (|3.4|). We also add a field 
independent term such that the action is finite at e — > 0. The pre-exponential factor D 
( |1.6| ) is a result of evaluating the path integral (|1 . 1| ) in the Gaussian approximation around 
the classical solution defined above, 



D 



me 



■TIE 



-2a 1 



Det' 



-d 2 + m 2 cosh(^)) 
—d 2 + m 2 



(3.6) 



where je = 0.577216 ... is the Euler constant. The first factor in ( |3.6j ) appears as a result 
of the mass renormalization. 
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3.1. Main semi-classical order 

We now proceed to the semi-classical calculation of the VEV ( e aLp ) using the repre- 
sentation ( p.!7| ). In order to apply the saddle-point machinery, it is convenient to begin 
with the form ( |2.6|) for the integral Xn- To the lowest semi-classical order, 

r coshf'y) 

Q[ 7 ] ~ e . (3.7) 

Here and below the notation 

r = mR (3.8) 
is used. The corresponding saddle-point equations have the form, 

r sinh(p fc ) = 2n (k + a) , k = 0, ±1, . . . , ±N . (3.9) 

In writing (|3.9|) we have assumed that b 2 — > and a = ab is fixed. Thus we find, 



1 N 

\ogZ N = -^ Yl cosh( Pfc )-2p fc (A; + a)}+0(l) . (3.10) 

k=-N 



The sum can be evaluated with the result, 

N 



{ - cosh(p fc ) -2p k (k + a) j = T - (—) log(— — j +2a 2 logiV+ 

(3.11) 

( 2iV(iV + 1) + i + 2a 2 ) log (^) + 4 ^ fc log(fc/e) - 4 log A G + O 



fc=-JV 

iV 

I 

fc=i 

where Ac = 1.282427... is the Glaisher constant and the function T = T(r, a) reads 
explicitly, 

/°° dr 
^ t sinh(r) 3?e [ log ( 1 - e~ r cosh M+ 2 — ) ] . (3. 12) 

-00 ^~ 

Now, combining Eqs. (|2.17| ), (|3.10| ), ( |3.11| ) one obtains 



(e a *) R ~e-l? , 

with 

S = S o (a)-T(r,a)+T(r,0) , (3.13) 

and 

2 1 ( m \ f°° dt f sinli (2at) 22/] , ^ 



Thus the function (|3.13| ) coincides with the regularized Euclidean action 



This result was obtained by a different method in 
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3.2. Semi- classical expansion of Q-function 

To compute the VEV to one-loop order, we have to find the next term in the semi- 
classical expansion of Q. It can be obtained by means of the TBA equation ( |2.13|) . The 
kernel $ ( p.!4j) allows an expansion in b 2 , 



m = 2n5{6) + 2nb* P.V. —jrf- + 0(6 4 ) , (3.15) 

sinn (6' J 

thus to lowest order the solution of the TBA equation has the form, 

e £ = e rcosh(e) - 1 + 0(b 2 ) . 

With this equation and the definition ( |2.15|) we calculate 

rcosh(7) rcosh(7) r7sinh(7) f°° dr log (l - e" rcosh(T) ) 2 



2ixb A 2 tv 2 7r J_ 00 2tv cosh(7 — r) 

(3.16) 

In order to see an analytical structure of this function it is instructive to represent Q in 
the form of an infinite product, 



1 , , < . r cosh(7) 
I r cosh(y) (rC ltj \ 2¥ , 

V2rx 



(3.17) 



n=l 

Here 7e is the Euler constant. 
3.3. One-loop order 

The saddle-point approximation allows one to find the one-loop order in the semi- 
classical expansion. To this order, 

N 

Tn = W J] sinh(p fc - P?n ) J] Q 2 [p k ] e 2 ^ a+ V (l + 0(6 2 )) . (3.18) 

N>k>m>-N k=-N 

Here the function W is a result of the Gaussian integrations in ( |2.6| ) around the saddle 
points 7 fc = p k (|J), 

yfc = -7V V rCOsh (Pfc) 
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The main steps in the calculation of ( |3.18| ) are given in Appendix B. Our final result has 
the form ( |1.6| ) with the function S given by ( |3.13| ) and 



log D = log D + T(r, a) - T(r, 0) - a d a T(r, a) - 

f d {r{ (9AT f - JL (a* rUo sit _ (8^r) 2 |„ ) } . (3 ' 19) 



Here 

, „ 2 , 1 f°° , sinh 2 (2crf) 

logD = -2« 2 log2 + - / dt ^ 

^ Jo cosh (t) 

Recall that ( ^.19| ) should coincide with the functional determinant 



4. High-temperature behavior 

Now that we have computed let us check the result for some limiting cases. Here 

we argue for R — > behavior of the VEVs. 

Due to the scaling properties of the interaction operator in one can rescale 

the problem to a circle of circumference 2tt. Thus, the Hamiltonian of the model under 
consideration takes the form 



H s hG = 



2n ^ 



J ^dr {A,lP + ±-{d rX T + » (§-) 2hq [J* + e-**) } , (4.1) 



where n = j j- is the momentum conjugate to \ = ip\ Xl= o- The mass of the sinh-Gordon 
particle is related to the parameter fx by [fL9l , 



n-b 2 ) 



7rr(l + 6 2 ) 



4^ 



(4.2) 



For r —>■ and a > the main contribution in the path integral ( |1 . 5| ) comes from a region 
of the configuration space corresponding to 

X ~ -2q log(r) > 1 . (4.3) 



In this region we can neglect the term e bx in the Hamiltonian ( |4.1| ), and approximate 
the ground state wave functional VPo[x] by a proper wave functional from the Liouville 
Conformal Field Theory (CFT). 
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More explicitly, the Hilbert space of the Liouville CFT contains a continuous set of 



primary states | p ) parameterized by p > with the conformal dimension ]20 

A 

We will assume that these states are canonically normalized 



2^ 



(4.4) 



(p' I V ) = 2tt 5(p — p') . 

Let ^ p [x] be a normalized wave functional corresponding to the state \p). As was discussed 
in |T7|,|18||, the following relation holds 



*oM ~L p y p [ X ] 



(4.5) 



in the region (^4.3| ). Here p = p(r) solves the equation 



p(r) : 2pq log 



8ttI (b 2 



bq 



71 i o 
h ksm 

2 



\og{T(l + 2ip/b)r(l + 2ipb) } . (4.6) 



We emphasize that L p is a unique coefficient provided a normalization of is chosen. 
Therefore, one can expect the following relation for r < 1: 



[e a r) R KL p L. 



(p\e aip \p) 



Liouv 



R (vac\vac) R 

The matrix element (p \ e a<p |p)Liouv was found in |2^,|18[. It reads explicitly, 



(4.7) 



R \ 2a(q — a) 

>k"*»Li.,»v= 



n/jT(b 2 ) b 



2\ u2-2b A 



r(i-6 2 ) 



-a/b 



T T(2a) T(2ip) ?(-2ip) 
T 2 (a)T(a + ip)T(a — ip) 



Here we use the notations [18 



logT(o) = -( £-a) log(26)+ / ^ 



q \ 2 -2t sinh 2 ( (qb — 2ab)t ) 
2~ a ) 6 ~ sinh(2t) sinh(2t6 2 ) 



and 



T = a a T(a)| a=Q . 
It is easy to see that the function p = p(r) (|4.6|) satisfies the condition, 



lim p(r) 
1 — >o 



(4.9) 
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Using Eqs.(|4.7| ) -( [19| ), we can derive 

-™r(^)r(i + ^i 2,,,, 



2q. 



where 



M 2 = 



^2aq 



R ,2a{a-q) T(2a)Tg 

T 4 (a) 



linip^o { 4p 2 L P L- P } 



(4.10) 



(4.11) 



R (vac\vac) R 

does not depend on a. In writing ( [4 . 1 0| ) we also used the relation (|4.2|). Unfortunately, the 
function M = A/"(r, b) is not known in closed form for an arbitrary value of the coupling 
constant. One can obtain its limiting value as b 2 — > 0. For b 2 <C r <C 1, it is sufficient to 
consider the dynamics of the zero-mode X [|TT|, 



2tt 



2^ 



X(r) 



(4.12) 



In this approximation, known as the mini-superspace approach 1 , the Hamiltonian (|4.1j ) 
is substituted by, 

H ms = |{-2d|+(^)\osh(&X)}. 
The Schrodinger equation 

H ms V (X) = E ms V (X) , 

coincides with the modified Mathieu equation and the wave functional \l/o is represented 
by its lowest eigenfunction. We will use the common normalization condition 



/oo 
dX ^l(X) = 1 
-oo 



(4.13) 



The Liouville wave functionals \& p in the mini-superspace approximation have the form, 

2 / r ^2i P /b 



F(2ip/b) \8nb 2 



K 2ip 



~ \4nb 2 



r bx 
e 2 



Here K u (z) is the MacDonald function. Now it is clear that the mini-superspace approxi- 
mation for M (|4.11| ) can be obtained from the large X behavior of the normalized Mathieu 
function V (X) (gig) , 



7T 



r cosh(6X/4) 



exp 



2-Kb 2 



cosh(6X/2) 



as 1^ ±oo . 



(4.14) 
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Of concern to us is the behavior of Af ms in the domain 6 2 C r C 1. In this case, we 
replace ^o(X) by its WKB asymptotic ( [4.14| ) and readily obtain, 

lim Af 2 = ^1 r < 1 . (4.15) 

b 2 ^0 {2-k) 6 

Having arrived at Eq. ( |4.15| ) , we can straightforwardly expand ( |4.10| ) in a power series 
of 6 2 , 

{e atp ) R nF (l + 0(6 2 )) (6 2 <r<l), (4.16) 



with 



T-> s 2 tt (a-l) 
/ It \ £3 



exp (4 5 (l/2 - a) - S (l/2 - 2a)) } : 



F= \2^J 2 ~~ — n.S (l (1./2-n)-,S (l CI./2-2n)) J-x 

o2 +^ — » „-7!Wr\§- 2 « /r(l-2a) T 2 (a) 



G V4tt^ y r(2a) T 2 (l-a) 
where Sq is given by Eq.(|3.14j) and Aq is the Glaisher constant. It is possible to show that 



the high-temperature expansion of (|1.6|) exactly matches Q4.16|) (see Appendix B for some 
details) . 

5. Low-temperature expansion 

We have mentioned in the Introduction that the finite volume VEVs can be understood 
as thermal averages ( |1.2|) . Hence, ( e a<p )r admits the low-temperature (R — > oo) expansion 
in the form, 

oo 

log((e«o) R /g a ) =l + J2G k (r) • (5.1) 

fe=i 

Here Gk represents /c-particle contributions in the infinite-volume channel and 

G k (r)~e- kr . 

Recently A. Leclair and G. Mussardo || proposed an integral representation which is 
sufficient to generate Gfc(r) systematically in terms of form-factors of the field e a<p at 
R = oo and the solution of the TBA equation ( J2.13| ). Taking into account contributions 
of one- and two-particle states to the trace (|1.2|) , they obtained 

/OO 
2~ f - {9) + 

J_ 00 2n 2-k cosh(6q - 6 2) 
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where the notations 

f-m = tt^) 

and 

sin 2 (^) 



are used. The function $ is the kernel in the TBA equation (|2.13| ). With (5.2) and the 



TBA equation one can calculate the first two terms in the low-temperature expansion ( |5.1| 



G, = 1M K (r) 

TV 



G 2 = 4[a] r ^i^i(^e 1 -e 2 )-2nS(e 1 -e 2 ))e- rcosh ^ +rcosh ^+ 

J-oo 27T 2tt V / 



(5.3) 



-00 

fCC 



o , d6^d02 ®(0l - O2) r cosh 0!+r cosh S 2 

L J ' 2tt 2tt cosh(6' 1 - 2 ) 



2 



where K n (r) is the MacDonald function. We now expand ( |5.3| ) as a power series in b 

G 1 = 4 j + s 2 (a) - a s(2a) + 0(b 2 ) j K (r) , (5.4) 

where 

sin(7TQ!) 

s(a) = . 

To expand G 2 one needs to use Eq.( |3.15p , 

G 2 ={ ^fi- + s 2 (2a) - 2a s(4a) \ K (2r)- 

l & J (5.5) 

4 s 2 (a) r 2 ( if?(r) - K 2 (r) ) - s 2 (2a) r 2 ( K 2 (r)K (r) - if 2 (r) ) + Q(b 2 ) . 



It is quite straightforward to verify that the low-temperature expansion of ([Oj) exactly 
reproduces (|5.4|) and 



6. Conclusion 

The proposed 7-representation ( |2.17| ) is the main result of this paper. Its rigorous 
derivation has not yet been achieved. Although (|2.17| ) are conjectures, the evidence pre- 
sented in this paper appears to make it reasonable to take them as the starting point for 
further investigation. 
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One can expect that similar representations exist for non-minimal CFT, say, the Li- 
ouville theory and SL(2,M)/U(1) non-compact a-model. It may cast new light on many 
unsolved problems of 2D Quantum Gravity. In this connection an intriguing similarity be- 
tween the integrals appeared in the 7-representation and Matrix Models of 2D Quantum 
Gravity |2^,£23],|25j can be mentioned. 
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Appendix A. 

The QFT defined by ( |1.3|) possesses infinitely many local IM I 2n -i whose vacuum 
eigenvalues have appeared in the equation ( |2.11| ). They can be represented in the form, 

hn-i = J 7^ ^T 2n (x 2 + ix u x 2 -ixi) + 6 2n -2(>2 + ix u x 2 - ix\) ) , 

where the local fields T 2n and 02n-2 satisfy the continuity equations, 

d B T 2n (z, z) = d z Q 2n _ 2 (z, z) . 

Although a general expression for the densities T 2ri , ©2n-2 is not known, they are deter- 
mined up to normalization by the commutativity conditions, 

[l2n-l j %.m-\ ] = . 



In Refs. HT5| , |27| the following normalization was adopted, 

T 2n = 2~ 2n (d z cp) 2n + . . . , (A.i; 
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where omitted terms contain higher derivatives of ip and exponential fields. Notice that the 
condition ( |A.1|) does not depend on the regularization scheme defining the composite field 
(d z (p) 2n . With normalization (|A.1| ) the constants C n (|2.11|) were found in Refc. [28 , [L5| , 



C r . 



r( 



(2n-l) b 
2q 



tv n\ q 



rn 



r(4)r(i) 



-I l-2n 



2qbi 



(A.2) 



The local IM l2n-i are certain deformations of the local IM of the Liouville CFT. Let 
^2n-i(f>) be an eigenvalue of the Liouville local IM corresponding to the state \p) 



while Ign-i is the sinh-Gordon ground state eigenvalues of l2 n -i- Eq.(O) suggests the 
following relation for r < 1: 



I 



2n-l 



<2n-l 



(p(r)) +0(r Aq \r Aq / h ) 



(A.3) 



where p(r) solves ( [L6| )li. The power corrections in r ( |A.3| ) can be obtained by means of 
Conformal Perturbation Theory. Explicit forms of the functions l2n-i(p) (f° r n = 1, . . . ,8) 
are given in Appendix B of Ref. |27]] . Here we present only the first two of them, 

T / N 27r ( 2 1 

I 3 (p) = 



2tt x - 



4 _ ^ 46 4 + 176 2 + 4 
P ~T + 96062 



In Tables 1-4 we list numerical values of the local IM l2 n -i for some 0.01 < r < 1 and 
b 2 = 0.81 which were obtained by means of numerical solution of the TBA equation ( |2.13j ) 
with use of Eqs.( |2.12"D and ( |A.2| ). These data are compared against values of the Liouville 
local IM Im-x (p( r ) )• We consider the content of Tables 1-4 to be an impressive evidence 
in support of the relation ( |2.12| ). 



For n = 1 this relation was discussed in Ref. | IS]. 
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r 


Ii 


h((p(r)) 


1.0 


-0.0059897196942029 


-0.0059891933248581 


0.8 


-0.0123637695005731 


-0.0123636397906571 


0.6 


-0.0183955019094376 


-0.0183954816200409 


0.4 


-0.0242191017573388 


-0.0242191003738075 


0.2 


-0.0301582132435655 


-0.0301582132312211 


0.1 


-0.0335250692109914 


-0.0335250692108919 


0.01 


-0.0381898149656469 


-0.0381898149656469 



r 


Is 


ia((p(r)) 


1.0 


0.01857760504756 


0.01858088002375 


0.8 


0.01975951264163 


0.01976027692024 


0.6 


0.02095100471765 


0.02095111804696 


0.4 


0.02216988492643 


0.02216989225147 


0.2 


0.02348269733500 


0.02348269739676 


0.1 


0.02425825249985 


0.02425825250033 



r 


Is 


h((p(r)) 


h 


J r ((p(r)) 


1.0 


-0.02830178173 


-0.02830149724 


0.0731035360 


0.0731032652 


0.8 


-0.02936366562 


-0.02936359281 


0.0750097817 


0.0750097140 


0.6 


-0.03040878256 


-0.03040877074 


0.0768644292 


0.0768644184 


0.4 


-0.03145618592 


-0.03145618508 


0.0787032289 


0.0787032281 


0.2 


-0.03256438579 


-0.03256438578 


0.0806285556 


0.0806285556 



r 


I 9 


h((p(r)) 


In 


hi((p(r)) 


1.0 


-0.29532264 


-0.29532199 


1.73235 


1.73235 


0.8 


-0.30116225 


-0.30116209 


1.75983 


1.75983 


0.6 


-0.30680503 


-0.30680500 


1.78627 


1.78627 


0.4 


-0.31236318 


-0.31236317 


1.81220 


1.81220 


0.2 


-0.31814537 


-0.31814537 


1.83907 


1.83907 



Tables 1-4. Comparison of the LHS and RHS of equation flAj) (b 2 = 0.81). 
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Appendix B. 



Here we proceed with calculation of the products in ( |3.18| ) and give some technical 
hints on the study of their high-temperature behavior. First, let us consider the product, 

e M * = W H sinh( Pfc - Pm ) . 

N>k>j>-N 



Using the relation 



d r Pk = ~- tanh(p fc ) 
r 



which follows from the saddle-point equations (|3.9|), one obtains 



d, 



1 N 

^ = 4 E 



cosh(/9 fc - p m ) 



2r k ^_ N cosh(p fc ) cosh(p r 



This sum can be rewritten in the form, 

(2/V + l) 2 I 



d r M 1 



2r 



+ 2r 



N 



^ tanh(p fc ) 



k=-N 



(B.l) 



Notice that the sum in ( P-l|) converges for N — > oo 



N 



lim tanh(pfc) = — -d r d a T + 2a . 



k=-N 



To derive this formula we used (|3.11| ) and the saddle-point equations ([T9|). Thus we obtain, 



Mi =M 



N 



(27V + 1) 



2cT lot 



2 

(^) -aS a T(r, «) - /°° y ^ (W) 2 + ^(iV- 1 ) . 



log 



47T 



(B.2) 



The constant -M;v here does not depend on r. To find how it depends on a, let us consider 
d a Mi. A similar calculation leads to the equation, 



1 2tt N 

d a M x = - ( rd r d a T - 4a ) ( d 2 T + 4 log(47riV/r) ) + — £ 



tanh(po) 



2tt 



4 

AT AT 

E E 

a=-iV fc=l 



/ <— ' cosh(po) + cosh(p m ) 

m= — N 



tanh (pfc(a)) 



tanh (p fc (-a)) 



cosh(p fc (a)) + cosh(p m ) cosh(p fc ( - a)) + cosh(p m ) 



;b.s) 
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It follows immediately from the last equation that 



d a Mi\ = -4a log - 



Therefore, we conclude that the constant M.^ in (|B.2| ) does not depend on a. Notice 
that Eq.( p.3| ) is very convenient for studying the high-temperature limit r —>■ 0. It is 
straightforward to show that for a > 0, 

/A 7 r\ 2N(N+1) o 



2 ^T{a),,^ I^Je 3 * 



n 



r(i+ a ) yr(i-a) fc AA 1 r(^ 1 -«) r (^ i + «) 

To finish the calculation of ( |3.18| ) one needs to evaluate the product, 



N 

,M 2 . 



e 



I I Q 2 [p fc ] e 2 ^ a+fc ) . (B.4) 



The first two terms of the semi-classical expansion for Q are given by ( |3.16| ). With the 
saddle-point equation (|3.9|) M 2 in ( |B.4|) can be written as, 



M 2 = M' 2 + M' 2 ' , (B.5) 

where 



N 

M^ = -(b~ 2 + l) cosh(p fc )-2p fc (A: + a)} (B.6) 



fc=-JV 

and 



J log(l _ e -^M T)) ^__. (B . 7) 
The sum M2 is evaluated by means of Eq. ( |3.11| ) . To calculate M'<1 we note that 

N 



t _ v cosh(6» fc ) y_ ( 



fc=-JV 

With the relations 

N - N 



2tv e rcosh ( T ) - 1 



E -4n = f E w* = f |a^ + 4iog(^ ))+o(i) 
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and 



one obtains, 



dr 



- dlT\ 



3h(r) _ I 4 ° la=0 



M» 



4nN 



16tt 2 q |a-u 4tt 2 
We specify the integration constant here using the condition, 



(B. 



M. 



2 |1 — >oo 







which follows from the definition (|B.7|) . The function T (|3.12|) satisfies Laplace's equation 



r 



- 1 d r (rd r T) + —dlT = 0. 



This allows one to calculate the second integral in ( |B.8| ). Thus we find, 

dr „ „, „o„ „, „ „, . /AtiN 



I J^ r d l T U=od 2 a T + T\ a=0 + rd r T\ a=0 log 



Finally we note that the most efficient way to study the r — > limit of Mi (P-4|) is based 
on the representation ( |3.f 7| ). It shows that for r — > with rcosh(7) fixed, 



r cosh (7) 
2^ 



(f + 0(6 2 )) 



V2fcosh(7/2) V47T- 
Hence, examination of the product (|B.4j) at this limit creates no difficulties at all. 
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